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Global properties of the Skyrme-force-induced nuclear symmetry energy 
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Large scale calculations are performed to establish the global mass dependence of the nuclear 
symmetry energy, a syrn (A), which in turn depends on two basic ingredients: the mean-level spacing, 
e(A), and the effective strength of the isovector mean-potential, k(A). Surprisingly, our results 
reveal that in modern parameterizations including SLy4, SkO, SkXc, and SkP these two basic 
ingredients of a sym are almost equal after rescaling them linearly by the isoscalar and the isovector 
effective masses, respectively. This result points toward a new fundamental property of the nuclear 
interaction that remains to be resolved. In addition, our analysis determines the ratio of the surface- 
to- volume contributions to a sym to be ~1.6, consistent with hydrodynamical estimates for the static 
dipole polarizability as well as the neutron-skin. 
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The knowledge of the nuclear equation of state (eos) 
for neutron-rich systems is of fundamental importance 
for nuclear physics and nuclear astrophysics. The sta- 
bility of neutron-rich nuclei, the r-process nucleosyn- 
thesis, the structure of neutron star, and the simula- 
tions of supernovae-collapse depend sensitively on the 
EOS and, in particular, on the nuclear symmetry energy 
(nse). Our study gives new insight into the two ba- 
sic ingredients of the symmetry energy, the mean-level 
spacing, e(A), and the effective strength of the isovec- 
tor mean-potential, n(A), revealing the striking property 

^■e(A) w ^k(A) that holds for velocity-dependent in- 
teractions having isoscalar and isovector m\ effective 
masses, respectively. This apparently fundamental prop- 
erty of the effective nuclear force remains to be explained. 

In asymmetric infinite nuclear matter (inm), in the 
vicinity of the saturation density, po, the EOS (the energy- 
density per particle) is conveniently parametrized using 
the following Taylor expansion: 



A 



K- 



-ay 



— (P-Po) 
Po 



18 Po 
AK, 



2 (p~ Pof + 



(1) 



ISA 



2-0- Po) 2 



where I = \N — Z\/A. The isoscalar inm saturation den- 
sity po, and the values of the volume binding energy ay, 
the incompressibility parameter Koo, and the asymme- 
try energy a sym serve as primary constraints for micro- 
scopic nuclear models. For modern Skyrme force param- 
eterizations which are subject of the present work pq w 
0.16fm" 3 , a v w -15.9 ± 0.2MeV, a sym « 32 ± 2MeV, 
and Koo ~ 225 ± 25 MeV. Higher-order curvature correc- 
tions to the NSE, p, AKoq, are rather poorly constrained. 
All these values are derived from the studies of finite nu- 
clei. 

Integrating out the r-dependence from the energy- 
density leads then to the semi-empirical mass formula 



(ld) for the energy per particle which is conventionally 
written as: 
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where as and a sym are coefficients defining contributions 
from the surface energy and the surface part of the sym- 
metry energy, respectively. There is at present no con- 
sensus concerning the magnitude, A, as well as origin of 
the term linear in which is often called the Wigner 
energy. 

Another controversy exists concerning the surface con- 
tribution to the NSE. The values of the surface-to- volume 
/ahQn quoted in_the literature vary 
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ratio r S / v 

strongly. For example, Danielewicz [l| estimates it to be 
2-0 < r s / v < 2.8, the mass formula of Ref. @ yields 
rs/v ~ 1-6 while the hydrodynamical-type models that 
include properly polarization of the isovector density pre- 
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2 [i| which, according to Ref. , is consistent 
within 10 -j- 20% with Skyrme-force (sf) calculations. 

The main objective of this work is to study the symme- 
try energy within the Skyrme-Hartree-Fock (shf) model. 
Notwithstanding, that a deeper understanding of the 
symmetry energy is crucial in order to reach a consen- 
sus on the existing variety of SF parameterizations, or to 
constrain the coupling constants of a more general local 
energy density functional (ledf). Our results point to- 
ward a deeper relation between the average level spacing 
and the strength of the mean isovector potential which 
has not been addressed hitherto. The calculations pre- 
sented below also allow us to determine the surface-to- 
volume ratio of the shf symmetry energy. 

In our previous letter |5j we have demonstrated that 
the SHF symmetry energy behaves rather unexpectedly 
according to the formula: 

E ( s S y " F) = \e{A, T Z )T* + \k{A, T z )T(T + 1) , (3) 
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where e(A,T z ) w e{A) and n(A,T z ) w re(A) are fairly 
independent on T 2 , at least for T 2 > 8, and denote the 
mean-level spacing at the Fermi energy in iso-symmetric 
nuclei and effective strength of the isovector mean- 
potential emerging within the SHF, respectively. More 
precisely, k is related to the isovector part of the SF in- 
duced ledf (s-ledf) H(r) = Y,t=o 1 ^*( r ) : 



H t = CI p\ + C? p Pt A Pt + ClfHTt + C; J J\ 



Definitions of all local densities and currents p, r, J as 
well as the explicit expressions for coupling constants 
Ct can be found in numerous references and we follow 
the notation used in Ref. [fj. Due to the isoscalar- 
density dependence of the SF , the coupling constants 
Cf[po] °f the S-ledf are functionals of po, giving rise 
to the isoscalar rearrangement mean-potential Uq = 
J2t=o l IT^Pt' Since our procedure of extracting e and 
k involves setting the C\ = 0, see Ref. |J, part of the Uq 
related to the Cf is formally treated as being related to 
the isovector degrees of freedom. Note that this separa- 
tion is is consistent with the way the symmetry energy 
constraint is superimposed on the SF . It does not af- 
fect the generality of our approach and similar analysis 
can be performed also within the much wider class of the 
Hohenberg-Kohn-Sham ledf theories. 

In the present work we focus on the global mass depen- 
dence of the SHF values of e(A) and k(A) and perform a 
systematic calculation covering all even-even nuclei hav- 
ing 20 < A < 128 from N = Z to almost the neutron 
drip line. Coulomb and pairing effects are disregarded 
i.e. the emphasis is on the strong interaction acting in 
the particle-hole channel. The calculations are performed 
for a set of different SF parameterizations as the SkPhj , 
SkXc 0, Sly4 0, SkO P3, SkM* [ll|, and SIII fTl, 
using the SHF code hfodd of Dobaczewski et al. |&ll3|. 

The procedure used to extract e(A, T z ) and k(A, T z ) 
follows exactly the one outlined in Ref. First, we 
set all the isovector coupling constants C\ = in the 
S-LEDF (0} and extract e(A,T z ) by comparing calcu- 
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AE^{A,T z ) = \s{A,T z )T\ 



(5) 



In the next step, we compute the total SHF binding energy 
Eshf(A,T z ) and compare: 



AE SHF (A, T z ) — AE^hf (A, T z ) = ^k(A,T z )T(T + 1), 

(6) 

in order to determine k(A,T z ). 

For each A and small T z , the values of e(A, T z ) oscillate 
quite rapidly. However, they clearly tend to stabilize for 
T z > 8 where e(A,T z ) « e(A). The values of k(A,T z ) 
appear to stabilize faster and k(A,T z ) r* k(A) essentially 
already for T z > 4. It should be mentioned that in the 
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FIG. 1: The isoscalar effective mass scaled values of 
e*{A) = —^-s(A) (stars) and k(A) (circles) calculated using 
the SHF method with SLy4 parametrization. Open symbols 
denote e*(A) and k(A) averaged over all the calculated nuclei. 
Filled symbols mark smoothed values of e*(A) and k(A) cal- 
culated using a restricted set of data. Vertical arrows indicate 
major shell gaps. Note the strong influence of shell structure 
on e*(A) and the smooth behavior of R(A). 



case of the SkO parameterization the formula JBJ docs 
work only approximately. For this force we observe an 
enhancement in the linear, ~T. This effect is, however, 
much weaker than the analo gous effect found recently 
within relativistic mean field |l4| where it restores the 
E S ym ~ T(T + 1) dependence of the total NSE. 

For further quantitative analysis of the mass depen- 
dence of the nse we use the mean values of e(A) and 
R(A). These averages over T z at fixed A are calculated 
using the following restricted set of nuclei: T z > 4 for 
A = 20; T z > 6 for A = 24; and T z > 8 for A > 28. 
By using a restricted set of nuclei we smooth out both 
e(A) and k(A) curves in order to diminish the possible 
influence of shell structure. 

The global mass dependence of the two components of 
the symmetry energy, s and R are presented in Fig. ^ 
Although representative for the SLy4 parametrization, 
the figure shows several features which appears to be in- 
dependent of the type of the SF parametrization. These 
universal features include: {%) strong dependence of e(A) 
on kinematics (shell effects); (ii) almost no dependence 
of k(A) on kinematics; (Hi) clear surface (~ A l /3 ) de- 
pendence reducing the dominant volume term (~ ^) in 
both e(A) and R{A). 

Indeed, the values of e(A) show characteristic kinks 
close to double-(semi)magic A-numbers. These kinks are 
magnified when all the calculated nuclei are used (no 
smoothing) to compute e(A), but without affecting qual- 
itatively the overall profile of the curve. On the other 
hand, R{A) is almost perfectly smooth with barely visible 
traces of shell structure. It confirms our earlier conclu- 
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sion [5j that the gross features of the Skyrme isovector 
mean potential can be almost perfectly quantified by a 
smooth curve parametrized by a small number of global 
parameters. 

In the analysis of a leptodermous expansion of e(A) 
and k(A) we consider volume (V) and surface (S) terms: 



The values of the isoscalar-effective-mass-scaled expan- 
sion coefficients e v , e* s as well as the values of the 
isovector-effective- mass-scaled expansion coefficients k v , 
Kg are collected in Tab.[I] First of all, let us observe that 
the calculated value of e v « lOOMeV corresponds to the 
pure Fermi gas estimate efg- This result can be under- 
stood based on the analytical expression for the Skyrme 
force nse coefficient in the limit of symmetric inm, a^L, 
provided that the standard textbook formula is rewritten 
in the following way: 
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where CI and Cf define the isovector part of the S- 
LEDF, see Eq. (J3J. Eq. (JSJ) clearly separates the con- 
tributions from the isovector and the isoscalar part and 
relates the latter to the single particle energies in INM, 

2 2 

£ p = 2WI + s (P' £ p) = 2^7* > with a self-energy term, 
E(p, e p ), that describes the interaction with the nuclear 
medium incorporated into the isoscalar effective mass. 
Hence, Eq. (JHJ further supports our interpretation of 
the NSE strength. 

The most striking result of our analysis is the near- 
equality of « k* occurring for all modern param- 
eterizations, see Tab. [I] and Fig. Indeed, e* differs 
from R* only for old parameterizations like the SIII and 
SkM* . This result confirms the rather loose claims often 
appearing in textbooks that "the kinetic energy [efg] 
and the isovector mean-potential contribute to the a sym 
in a similar way" is indeed correct but only after disre- 
garding non-local effects. To our knowledge, it has never 
been discussed why this apparently independent quanti- 
ties should be similar. 

The relation e* w k* could be accidental. However, 
given the fact that these modern Skyrme forces have been 
fitted in a rather different manner, suggests that the near 
equality is of fundamental nature. To find an explanation 
of this result will be a challenge for further studies. Note 
that the analytical relation © does not show any explicit 
scaling due to the presence of the isovector effective mass, 
m* . This is probably not very surprising since m\ defines 
the enhancement in the energy weighted sum rule for the 
translational symmetry violating (finite nucleus) dipole 
mode. Moreover, the key isovector coupling constants 
C\ and Cf vary completely erratically from force to force, 
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FIG. 2: The correlation between the effective mass scaled 
volume Ey and kJ> (open dots) and the surface e$ and Kg 
expansion coefficients. Note, that except for the SIII and 
SkM* interactions, the expansion coefficients are equal, see 
also Tab.El 



although we have found that they are linearly correlated 
for all the SF studied here. It indicates that they are 
neither well understood nor well constrained. 

The advantage of our interpretation of e(A) becomes 
evident when evaluating r £ using the semi-classical ap- 
proximation |15| in straight forward fashion. The appro- 
priate formula which takes into account the diffuseness 
of the potential, see Ref. |Ty| : 
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where is the level density at the Fermi energy, kp w 

1.36 fm -1 while Vm and Sm denote volume and surface 
matter-distribution, respectively. Assuming spherical ge- 
ometry ^M- Rj - and adopting for r w 1.14 fm, i.e. 
the value consistent with the standard Skyrme force sat- 
uration density po ~ O.lGfm" 1 , one obtains r e « ~ 
1.52 which is indeed very close to the calculated ratios, 
see Tab. HI 

The SHF models yield r$/v ~ 1-6 in accordance with 
the LD ratio |2|. The static dipole polarizability (SDP) 
an [o~(lu) denotes photo-absorption cross-section]: 



f (T { UJ ) , r, 
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provides an independent cross-check of the ratio r S / V . 
Indeed, using the so called hydrodynamical model simple 
estimate for an can be derived [3j: 
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is the so called Migdal SDP 



value |17| which is valid for large systems with negli- 
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SkP 


1.000 


0.741 


108.8 


175.1 


1.61 


106.0 


163.1 


1.54 


1.03 


1.07 


0.95 


30.0 


31.5 


49.4 


1.57 


SkO 


0.896 


0.852 


107.2 


166.2 


1.55 


110.6 


176.1 


1.59 


0.97 


0.94 


0.79 


32.0 


31.2 


49.0 


1.57 


SkM* 


0.789 


0.653 


106.3 


180.9 


1.70 


71.4 


107.3 


1.50 


1.49 


1.69 


1.37 


30.0 


30.5 


49.2 


1.61 


SIII 


0.763 


0.655 


97.5 


143.8 


1.47 


75.2 


103.2 


1.37 


1.30 


1.39 


1.34 


28.2 


30.3 


43.3 


1.43 



TABLE I: The table includes: the isoscalar, mj/m, and the isovector, m*/m, effective masses; the volume, e v {k v ), the surface 
£ s ( K s)y an d the ratios r e = e* s /e v (r K = n s /n v ) of the expansion coefficients J7J of the isoscalar effective mass scaled s*(A); 
and the isovector effective mass scaled k* (A) = ^-R(A), respectively; the ratios of volume r v = e v /n v and surface r s = s s /k s 
expansion coefficients and the inm estimate r*^ = ef^/^ftx,); the inm estimate ai^m and the calculated values ai^m, aiyln 

and rs/v = asym/ai^L of the symmetry energy coefficient as defined in @. The values of £ V {S)j k V(S)' a syrn, aij/m, and aiym 
are given in MeV. 



gible surface contribution. Using (r 2 ) = |i? 2 A where 

R = l^A^Hm and a%} n = 30MeV one obtains, in 
the Migdal's limit, the following numerical estimate 
cr (M)^_5 /3 _ i.73^b/MeV. Using this estimate and the 

experimental value of a^A' 5 / 3 2.7 ± 0.2/ib/MeV 
which is almost constant for A > 100 0,0], one obtains 
rs/v ~ 1-65 what is again consistent with our estimates 
based on the SHF calculations. 

The neutron skin thickness is another quantity which 
sensitively depends on isovector properties. Neutron rms 
radii are even today rather poorly known and can there- 
fore not be used to constrain the ledf or the EOS, see 
however 0] . For heavy nuclei one can continue to apply 
the hydrodynamical model of Ref. Q and evaluate the 
neutron skin thickness: 

*f.„^/i+ar^i- 1 (12) 

<r 2 > A \ i+ 3AV3 [iZ) 

Interestingly, Eq. i|12|) does not depend on the 
bulk NSE coefficient but only on the ratio rg/ v . For- 
mula l|12f> taken at rs/v ~ 1-65 gives again a very con- 



sistent results with our microscopic SHF calculations. 

In summary, the global mass dependence of 
the NSE strength a sym (A) and its two basic ingredi- 
ents related to the mean-level spacing, e(A), and to the 
mean-isovector potential, k(A) is studied in detail within 
the SHF theory. Our interpretation of the symmetry en- 
ergy enables us to unambiguously establish the surface- 
to-volume ratio of a aym (A), rs/v ~ 1-6 m agreement 
with the LD value of Ref. (J- This ratio is consistent 
with simple hydrodynamical estimates for the SDP and 
neutron skin thickness. The most striking results of our 
calculations is the near-equality of e* rj k* revealing that 
contribution to a sym due to the mean-level spacing and 
due to the mean-isovector potential are similar but only 
after disregarding non-local effects. Whether this is a 
fundamental property of the nuclear mean field is an open 
question that requires further studies. 
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